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Abstract 

We calculate dilepton production rates from a fireball adapted to the kinematical conditions 
realized in ultrarelativistic heavy ion collisions over a broad range of beam energies. The freeze-out 
state of the fireball is fixed by hadronic observables. We use this information combined with the initial 
geometry of the collision region to follow the space-time evolution of the fireball. Assuming entropy 
conservation, its bulk thermodynamic properties can then be uniquely obtained once the equation 
of state (EoS) is specified. The high-temperature (QGP) phase is modelled by a non-perturbative 
quasiparticle model that incorporates a phenomenological confinement description, adapted to lattice 
QCD results. For the hadronic phase, we interpolate the EoS into the region where a resonance gas 
approach seems applicable, keeping track of a possible overpopulation of the pion phase space. In this 
way, the fireball evolution is specified without reference to dilepton data, thus eliminating it as an 
adjustable parameter in the rate calculations. Dilepton emission in the QGP phase is then calculated 
within the quasiparticle model. In the hadronic phase, both temperature and finite baryon density 
effects on the photon spectral function are incorporated. Existing dilepton data from CERES at 158 
and 40 AGeV Pb-Au collisions are well described, and a prediction for the PHENIX setup at RHIC 
for Vs = 200 AGeV is given. 



1 Introduction 

There are convincing arguments from lattice simulations jlj] that QCD undergoes a phase transition at 
a temperature of about 150 - 170 MeV || from a confined hadronic phase to a phase where quarks and 
gluons constitute the relevant degrees of freedom, the quark-gluon plasma (QGP). It is believed that 
chiral symmetry restoration, associated with the vanishing of the chiral condensate (qq), takes place at 
about the same temperature. This new state of matter is clearly of enormous interest, and the search 
for experimental evidence of its creation in ultrarelativistic heavy ion collisions (URHIC) is going on. 
However, it is still a matter of debate whether the QGP has already been encountered. Any signature of 
the QGP is folded with the time evolution (and consequently the evolution of volume, temperature and 
baryon density) of the fireball created in such a collision. Furthermore, this evolution continues after 
the phase transition and any information is thus mixed with signals from the hadronic phase. Therefore, 
extracting an unambiguous signal for the QGP phase is not an easy task. 

Dileptons (e + e~ and fi + fj," pairs) are interesting probes in this context since they do not interact strongly 
and escape unthermalized from the hot region at all stages of the collision. Therefore, in contrast to 
hadronic probes which can tell only of the late stages (the freeze-out), dileptons carry information also 
on the early moments of the collision where the QGP phase is expected to exist. In the QGP phase, 
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dilcptons originate mainly from qq annihilation processes, whereas in the hadronic phase the main sources 
are pion and kaon annihilation processes which are enhanced due to the formation of the light vector 
mesons p, lo and <f>. In order to interpret signals from the QGP, one must obviously first understand the 
hadronic phase. 

Unfortunately, the initial conditions of a fireball formed in URHIC are not well under control, nor is 
its subsequent evolution known in sufficient detail so far. Different scenarios (pure hadronic, pure QGP 
and mixed phases) can all be designed such that they are in rough agreement with experimental data. 
This is shown in calculations using microscopic transport equations (j, P , hydrodynamical descriptions 
|^|, 0] and simplified approaches g, g|. In this paper, we make attempts to tighten the constraints on the 
QGP equation of state (EoS) and the fireball evolution dynamics by establishing contacts with lattice 
QCD thermodynamics and hadronic observables. We extend the work of |n| (and other previous work) 
in several directions. 

First, we develop a thermodynamically consistent quasiparticle picture, incorporating important aspects 
of confinement, to describe both the fireball dynamics and the dilepton radiation from the quark-gluon 
phase at the early stages of the collision. Perturbation theory is known not to work in the quark-gluon 
phase (l| even at temperatures well above Tq. We adopt instead an approach which treats quarks and 
gluons as massive thermal quasiparticles, with their properties determined to reproduce lattice QCD 
results. The driving force of the phase transition, the confinement /deconfinement process, is given a 
statistical meaning in terms of a suppression of thermally active partonic degrees of freedom as the 
critical temperature is approached from above. For the hadronic phase, we calculate the dependence of 
the photon spectral function (that enters the dilepton rate) on temperature and baryon density using 
Vector Meson Dominance combined with Chiral Dynamics. 

Second, a fireball model reminiscent of hydrodynamics is set up. The endpoint of the fireball evolution 
is fixed by requiring agreement with hadronic observables, such as Hanbury-Brown-Twiss (HBT) radii, 
mt-spectra and dN/dy distributions, which carry information on the freeze-out stage. Information on the 
initial state can be gained by estimates of the nuclear overlap geometry. We then construct the evolution 
between initial and freeze-out conditions by following the fireball thermodynamics along volume slices 
of constant proper time r. This is done under the assumption of entropy conservation where the total 
entropy of the system is determined at freeze-out. The EoS in the QGP phase is calculated using 
a quasiparticle picture. In the hadronic phase a smooth interpolation to a hadronic resonance gas 
approach is applied. We thus obtain the temperature T of the system at a given proper time. This 
procedure results in a thermodynamically self-consistent description of the fireball evolution, which is 
fixed by hadronic observables only. 

Using this model as an input, we calculate dilepton rates and compare results to existing CERES/NA45 
data of 40 and 160 AGeV Pb-Au collisions, finding good agreement in both cases. We also give predictions 
for the kinematical conditions at RHIC. Differences between the present approach and other works are 
discussed. 



2 Dileptons from a fireball 

The lepton pair emission rate from a hot domain populated by particles in thermal equilibrium at 
temperature T is proportional to the imaginary part of the spin-averaged, photon self-energy, with these 
particles as intermediate states. The thermally excited particles annihilate to yield a time-like virtual 
photon with four-momentum q which decays subsequently into a lepton-antilepton pair. The differential 
pair production rate is given by 
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where a — e 2 /4ir, j3 — 1/T, and we have neglected the lepton masses. We have defined H(q) = — IP /3 
and introduce the averaged photon spectral function R(q, T) — (12ir/q 2 ) Imfl(g, T). Here II^ denotes the 
trace over the thermal photon self-energy which is equivalent to the thermal current-current correlation 
function 

lWg,T) = ijd A x e^(T3^x) Ju (0))p, (2) 

where is the electromagnetic current. Eq.([j]) is valid to order a in the electromagnetic interaction 
and to all orders in the strong interaction. 

In order to compare with experimental data, we construct a model for the space-time evolution of a 
heavy-ion collision, assuming approximate thermal equilibrium to be a useful concept 113] . Some recent 
discussions suggest jlj], [l5| that equilibration times at the conditions of heavy ion collisions at SPS and 
RHIC may indeed be very short (possibly less than 1 fm/c), small compared to expansion times of order 
10 - 20 fm/c, although this is still a matter of debate. The proper approach for our purposes is to use a 
simplified fireball model which employs a hot cylinder, expanding isentropically both in the longitudinal 
and transverse direction. Temperature and density are assumed to be spatially homogeneous at a given 
proper time. In chapter ^[ this fireball model is explained in detail. 

The differential rate of eq.Q is integrated over the space-time history of the collision to compare the 
calculated dilepton rates with the CERES/NA45 data Q taken in Pb-Au collisions at 160 AGeV 
(corresponding to a cm. energy of */s ~ 17 AGeV) and 40 AGeV (y/s ~ 8 AGeV). The CERES 
experiment is a fixed-target experiment. In the lab frame, the CERES detector covers the limited 
rapidity interval rj = 2.1 — 2.65, i.e. Ary = 0.55. We integrate the calculated rates over the transverse 
momentum px and average over r/, given that d A p = Mpx dM dr\ dpr d0. The formula for the space-time- 
and p-integrated dilepton rates hence becomes 
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where tj is the freeze-out proper time of the collision, V(t),T(t)) describes the proper time evolution of 
volume elements moving at different rapidities and the function Acc(M,t],pt) accounts for the experi- 
mental acceptance cuts specific to the detector. At the CERES experiment, each electron/positron track 
is required to have a transverse momentum pt > 0.2 GeV, to fall into the rapidity interval 2.1 < r\ < 2.65 
in the lab frame and to have a pair opening angle O ee > 35 mrad. Eq.(||) is then convoluted with the 
finite energy resolution of the detector. Finally, for comparison with the CERES data, the resulting rate 
is divided by dN c h/dr), the rapidity density of charged particles. 

RHIC operates as a collider experiment, so in this case the fireball is centered around rj = 0. Here, the 
PHENIX detector acceptance can be schematically modelled by requiring that each electron/positron 
track falls in the rapidity interval —0.35 < rj < 0.35, has transverse momentum p? > 0.2 GeV and a pair 
opening angle of e e > 35 mrad. At present, an abundance of data on Au-Au collisions at = 130 
AGeV have already been analysed, and first data of the run at the higher energy y/s — 200 AGeV are 
available. 



3 Quasiparticle description of the quark-gluon phase 

In order to proceed, we need to develop a picture of the QCD dynamics in the quark-gluon plasma phase 
close to the phase transition where the system presumably spends an appreciable fraction of its time in 
URHIC at current kinematical conditions. A realistic equation of state is needed to model the expanding 
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fireball, and the calculation of the differential dilepton emission rate requires the photon self-energy with 
partonic intermediate states. 

Since it is expected that a s (T) — > at very high temperatures, one might naively expect that a perturba- 
tive approach to calculate these quantities is possible. However, it is now well-known that perturbation 
theory breaks down for gauge theories at finite temperature because of the occurrence of infrared diver- 
gences [|ll| . The Hard Thermal Loop (HTL) calculation scheme consistently resums all contributions to 
a given order in the gauge coupling constant g s and is explicitly gauge invariant. Its application relies 
on a scale separation into 'hard' (p ~ T) and 'soft' (p ~ g s T) momentum modes with g s <C 1, though 
for QCD, this condition is likely to be fulfilled only at extremely large T far outside the scope of present 
and future experiments fL2| . In addition, non-abelian gauge theories also suffer from a dependence on 
an 'ultra-soft' magnetic mass rn mag ~ g 2 s T which is intrinsically nonperturbative. 

The imaginary part of the photon self-energy Imfi appearing in eq.(|l|) has been calculated to two- loop 
order in thermal perturbation theory jl7j , and there exist estimates for the three-loop contributions Jul . 
So far, the convergence of the resulting series is badly behaved even for small couplings g s . A similar 
behaviour has been found in the perturbative calculation of the free energy of a QGP and also in simple 
scalar theory. Furthermore, close to Tq we expect intrinsically non-perturbative dynamics to enter: the 
confinement/deconfinement transition and chiral symmetry breaking are not amenable in an expansion 
in g s . In the light of these facts, we will use a more phenomenological approach to QCD thermodynamics 
in this work. 

One way to obtain nonperturbative input is through finite temperature lattice simulations of a pure 
gluonic or quark-gluon plasma. We have shown recently that it is possible to describe the equation of 
state (EoS) of such systems to a very good approximation by the EoS of a gas of quasiparticles with 
thermally generated masses, incorporating confinement phenomenologically by a temperature-dependent 
effective number of active degrees of freedom. Here we give a short summary of the method and refer 
the reader to |L9| for a more detailed discussion. 

From asymptotic freedom, we expect that at extremely high temperatures the plasma consists of quasifree 
quarks and gluons. HTL perturbative calculations find, for thermal momenta, spectral functions of the 
form S(E 2 ~k 2 — m 2 (T)) with m(T) ~ g s T. As long as the spectral functions of the thermal excitations at 
lower temperatures resemble qualitatively this asymptotic form, a quasiparticle description is expected 
to be applicable. The QCD dynamics is then incorporated in the thermal masses of the 'effective' quarks 
and gluons and in a function B(T) which is necessary for thermodynamical consistency and plays the 
role of the thermal vacuum energy. 

Since the existence of a preferred frame of reference (in which the hot matter is at rest) breaks Lorentz 
invariance, new partonic excitations (longitudinal gluonic plasmons and helicity-flipped plasminos) are 
also present in the plasma. However, their spectral strengths are exponentially suppressed for high tem- 
peratures and momenta k ~ T which dominate the integrals of macroscopic thermodynamic quantities 
such as pressure, entropy and energy density, therefore their contributions can be neglected. This is also 
in line with results from a lattice calculation of the quark and gluon propagators |l6| . For the differential 
dilepton rate, their influence is not so easily discarded, and we will discuss that in section || 

The thermal excitations can then be described by a dispersion equation 

E 2 (k,T) = k 2 +m 2 {T). (4) 

Here, k = |k|, and the subscript i labels the particle species: i = g for gluons and i = q for quarks. 
mj(T) stands for a thermal mass which is derived from the self-energy of the corresponding particle, 
evaluated at thermal momenta E,k ~ T. This is expected to be meaningful as long as the self-energy 
is only weakly momentum dependent in that kinematic region. Additionally, for a quasiparticle to be a 
meaningful concept at all, we require the imaginary part of the self-energy, and hence its thermal width, 
to be small. 
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The gluon mass, following ref. becomes 



T 



with the effective coupling specified as 



g(T,N c ,N f ) 



9o 



2N 
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j T 



(5) 



(6) 



Nc and Nf stand for the number of colours and flavours, respectively. The functional dependence of 
m g (T) on T is based on the conjecture that the phase transition is weakly first order or second order 
which suggests an almost powerlike behaviour m ~ (T — To) 1 with some pseudo-critical exponent 7 > 0. 
Setting go = 9.4, S = 10~ 6 and 7 = 0.1, the effective mass (||) approaches the HTL result at high 
temperatures. 

The thermal quark mass reads 



m qi T ) 

T 



\ 




N 2 - 1 



(V 



with the zero-temperature bare quark mass to 9i o- Note that in contrast to previous quasiparticle models 
extrapolated from HTL calculations [H], the thermal masses used here drop as Tc is approached 
from above. Of course, for T ^S> Tc, the near-critical behaviour inferred in (|^) ceases to be valid and the 
perturbative limit of m g (T) and m q (T) will be recovered. 
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Figure 1: The confinement factor C(T) as a function of temperature for two light flavours (m u ,d = 0) 
and one heavier flavour (m s ~ 170 MeV). 
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The pressure of the QGP system takes the following form: 

14 N f n, r 

p(T) = -JL / dk [C(T)f B {ED) dk [C(T)f D (ED] — - B(T). (8) 

Energy density e(T) and entropy density s(T) follow accordingly (see Oil for details): 



f°° . . oat r°° 

e(T) = -^/ dkk*[C{T)f B {El))El+Y,^T dkk 2 [C(T)f D (El)]El + B(T), (9) 
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s{T) = Ykl dke[C{T)f B {Ei r k +7^ +^^ f dkP[C(T)f D (El)] 



(10) 



Here, i?^ = k 2 + m 2 g {T) and = y/k 2 + m 2 (T) for each quark flavour q = i. The multiplicity 

^ s = f 6 counts the transverse gluonic degrees of freedom. The bosonic particle distribution function 
reads f B (E) = (exp(E/T) - l)" 1 , the fermionic one is f D (E) = (exp(E/T) + l)" 1 . 

The function C(T) accounts, in a statistical sense, for the onset of confinement by reducing the number 
of thermally active degrees of freedom. Its explicit form is obtained by calculating the entropy density 
of the QGP with the gluon mass (||) and the quark mass (0). Dividing the lattice entropy density by 
this result yields C(T). It can be parametrized as 



C(T) = C ([l + <y-^) • (Hi 



7c 

T 

For two light quarks and one heavy quark, the parameters take the values Co = 1.16, 5 C = 0.02 and 
7 C = 0.29 (see figure 0)0. The function B(T) is now uniquely determined from m^T) and C(T) up to 
an integration constant that is fixed by Gibbs' condition pqgp = Phadr a t Tq. 

Figure ^| shows the pressure, energy and entropy density for two light quark flavours (m u .d = 0) and a 
heavier strange quark (m s ~ 170 MeV) in our confinement model. 

We now proceed as follows: We adopt the picture of the QGP as a quasiparticle gas consisting of massive 
'quarks' and 'gluons' which are the relevant thermodynamical degrees of freedom. The resulting EoS 
serves as input for the construction of the fireball which is discussed in section [j| Furthermore, within this 
framework the quasiparticle qq-loop becomes the only contribution to the photon self-energy in eq.(jl|) 
since the quasiparticles are by construction non-interacting. The reduced probability for thermally 
exciting a quark because of confinement, and the thermal effective masses, have to be taken into account 
properly as discussed in the next section. 



4 Calculation of the photon spectral function 



t^We mention that the proposed method relies on input from the lattice. Whereas high precision data exist in the 
pure gauge sector, calculations with dynamical quarks are not yet in a satisfactory position to yield proper continuum- 
extrapolated results with physical quark masses. We estimate that the results obtained within our confinement model may 
still change by 5-15% in the vicinity of Tq once high statistics data are available. However, this small correction does not 
influence the results of the following discussion. 
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Figure 2: Pressure, energy density and entropy density for two light quark flavours (m UjC 2 = 0) and a 
heavier strange quark (m s ~ 170 MeV) in our quasiparticle model. The arrow indicates the ideal gas 
limit of massless three-flavour QCD. The function B(T) is also shown. 



4.1 The quark-gluon phase 

As long as the thermodynamically active degrees of freedom are quarks and gluons, the timelikc photon 
couples to the continuum of thermally excited qq states and subsequently converts into a charged lepton 
pair. The calculation of the photon spectral function at the one-loop level is performed using standard 
thermal field theory methods. The well-known leading-order result for bare quarks and gluons as degrees 
of freedom is: 




where q = (q°, q) is the four-momentum of the virtual photon, e/ the quark electric charge and m/ the 
quark mass of flavour /. This result, however, holds only up to perturbative higher order corrections 
in g s that take into account collective plasma effects. Here, contributions from soft gluons lead to 
strong modifications. The corresponding two- and three-loop contributions show no clear convergence 
Jl7| . [f8| . Close to the phase transition, we also expect non-perturbative confinement physics to enter. 
Consequently, we follow a different approach. 

Recalling the results of the previous section, the thermodynamic properties of the QGP as given by 
lattice QCD are well reproduced by a gas of quasiparticles. Let us now assume that a quark quasipar- 
ticle couples to a photon in the same way as a bare quark (a form factor representing the 'cloud' of 
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the quasiparticle could in principle also be included, but in absence of information about the detailed 
quasiparticle structure we ignore this point). For a gas of non-interacting quasiparticles, the one-loop 
result for Imll is adequate, with input properly adjusted. All higher order QCD effects manifest in the 
thermal quasiparticle masses m(T), the function B(T) and the confinement factor C(T). Incorporation 
of the first two features in the calculation is straightforward. The bare quark masses in eq. ( JT^ ) simply 
have to be replaced by the T-dependent quasiparticle masses for each flavour, see eq. fjj). The thermal 
vacuum energy B(T) does not contribute to the dilepton rate. 

The naive replacement fn — * C{T)fu is, however, not permitted in eq. (|l^). Since any modification of 
the free particle distribution functions leads to non-equilibrium field theory, products of delta functions 
(pinch singularities) may arise in loop calculations. Therefore, the quasiparticle model as it stands 
cannot be used in expressions derived from simple perturbative thermal field theory. Recalling the 
physical interpretation of the confinement factor C(T), we can use the expression for the dilepton rate, 
eq. (|l|) , instead. The mechanism for dilepton production at tree- level is the annihilation of a qq pair into 
a virtual photon where the quark lines are multiplied by the distributions /d(T), giving the probability 
of finding a quark or an antiquark in the hot medium. This also becomes clear when we look at the 
limit q -> of eq. ©. Then, 

lmU(q°,T) ~ Imn(q°,T= 0) • (1 - 2f D (q°/2)), (13) 

and the temperature enters only in the Pauli-blocking of the quarks propagating in the loop. Now, from 
eq. (0) 

dN 



f B (q u )ImU(q,T). (14) 



d 4 xd 4 q 

Combining the different thermal occupation factors, we end up with the well-known result 

^~[/^°/2)] 2 , (15) 

so the differential dilepton rate is proportional to the probability of finding a quark q times the proba- 
bility of finding an antiquark q with the correct momentum, as anticipated^]. The incorporation of the 
confinement factor is now obvious: since it reduces the number of thermally active degrees of freedom, 
it also reduces the dilepton rate by a factor of C(T) 2 . 

In summary, eq. (|]) can be used to calculate the dilepton rate originating from a hot QGP phase, 
provided an overall factor C(T) 2 is applied to account for the reduced probabilities, and the bare masses 
rrif in the one-loop expression (^) are replaced by the T-dependent thermal masses (Q). The role of the 
factor C(T) is illustrated in figure |^, where the differential dilepton rate originating from a hot QGP 
in the quasiparticle approach is shown for different temperatures. Note that the plotted quantity is 
independent of the fireball volume, so the resulting differences are only due to the dropping quasiparticle 
masses and the squared confinement factor C(T), which is responsible for a decrease by more than an 
order of magnitude at T ~ Tq as compared to the highest temperature shown. One also observes that, as 
expected, the (negative) slope of the production rate in the region of high invariant mass gets steeper as 
the temperature decreases. It is important to note that this setup neglects contributions from hadronic 
degrees of freedom above Tq- As mentioned, quarks and gluons become clustered into hadrons (glueballs, 
mesons) as the temperature approaches Tq from above. These hadronic excitations are comparatively 
heavy and thus do not contribute much to the thermodynamics. Since we do not know in detail how the 
statistical re-arrangement of degrees of freedom occurs, we refrain from including these hadronic sources 
of dilepton yield above Tq. Our calculation is therefore expected to give a lower limit on the leptonic 
radiation from the QGP phase. 

t'We neglect a possible chemical potential for the quarks. For a finite ^t, the corresponding expression would be 
3^ ~ - ^)/ 2 )M(9° + M)/2). 
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Figure 3: Dilcpton rates originating from the QGP phase in the quasiparticle model for different tem- 
peratures. 



The quasiparticle model does not take into account the collective plasma modes, such as the (longitudi- 
nal) gluonic plasmons and the (hclicity-flipped) quark plasminos. Since their residues arc exponentially 
suppressed in the HTL approximation for thermal momenta k ~ T ', their contributions to the thermo- 
dynamical quantities are negligible. However, in the case of soft dilepton production it is well known 
that these modes lead to sharp, distinct structures in the spectrum, referred to as Van Hove singularities 
p3| . The plasmino branch has a minimum in its dispersion relation at k 7^ (which follows on very 
general grounds |24j , independent of the HTL approximation) . This leads to a diverging density of states 
which, in turn, shows up in the dilepton spectrum as a pronounced peak. Our model cannot exhibit, by 
construction, such plasmino effects. However, since the peaks are roughly located at ~ 2m q (T), where 
m q (T) is the thermal mass of the quasiparticles, these Van Hove singularities would be smeared out by 
the fireball evolution. The thermal mass drops as the temperature goes down, dragging along the peak 
position with it. Furthermore, since m q (T) is of the order of the temperature T or smaller in our model, 
the singularities appear at low invariant mass (< 500 MeV) where they are overwhelmed by the hadronic 
part of the dilepton production. Therefore the presence of these collective modes would presumably not 
influence our results for the dilepton rate in the QGP phase. 

Our model of the QGP phase also finds preliminary support from a recent lattice calculation |5lJ of 
thermal vector meson correlation functions above Tq in quenched QCD with Clover improved Wilson 
fermions. Using the maximum entropy method, the vector spectral function was extracted from the 
corresponding current correlator. Although the statistical uncertainties are still considerable, it is inter- 
esting to note that the resulting spectrum resembles the free spectral function, as in our case, and has a 
gap at low energies given by a thermal mass threshold of (2 — 3)T, which is indeed close to 2m q (T), the 
natural cut-off of the spectrum and, correspondingly, of the thermal dilepton radiation in the confine- 
ment model (see figure [}]). Furthermore, this result seems to rule out heavy quark quasiparticles in the 
deconfined phase, as predicted by other phenomenological models J2^, ^l). Of course, higher statistics 
and improved actions are mandatory to confirm these observations. 
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4.2 The hadronic phase 



Below Tc, confinement sets in and the effective degrees of freedom change to colour singlet, bound qq 
or qqq (qqq) states. The photon couples now to the lowest-lying 'dipole' excitations of the vacuum, the 
hadronic J p = 1~ states: the p, u and <j> mesons and multi-pion states carrying the same quantum 
numbers. The electromagnetic current-current correlation function can be connected to the currents 
generated by these mesons using an effective Lagrangian which approximates the SU(3) flavour sector of 
QCD at low energies. The appropriate model for our purposes is the improved Vector Meson Dominance 
model combined with chiral dynamics of pions and kaons as described in p5| . Within this model, the 
following relation between the imaginary part of the irreducible photon self-energy ImH and the vector 
meson self-energies Tlv(q) in vacuum is derived: 

9 \ „2 



Irnn(,) = £^^|iV( 9 )| 2 , F v (q) = a ^ f[ T J""" (16) 
y 9v q 2 -m v + ilmll v (q) 

where my are the (renormalized) vector meson masses, g v is the 7V coupling and g is the vector meson 
coupling to the pseudoscalar Goldstone bosons n^, ir° and K ± , K°. Eq.([j"6|) is valid for a virtual photon 
with vanishing three-momentum q. For finite three- momenta there exist two scalar functions H^ and 
ILt, because the existence of a preferred frame of reference (the heat bath) breaks Lorentz invariance, 
and one has to properly average over them. However, taking the limit |q| — > should be reasonable for 
our purposes in view of the fact that the cm. rapidity interval accessible at CERES and RHIC restricts 
|q| on average to only a fraction of the vector meson mass my. 

Finite temperature modifications of the vector meson self-energies appearing in eq.(|l6|) are calculated 
using thermal Feynman rules. The explicit calculations for the p- and </>-meson can be found in ref . [[To| . 
At the one-loop level, the p and 4> are only marginally affected by temperature even close to Tc because 
of the comparably large pion and kaon masses: ~ Tc, mj< — 3 Tc- The thermal spectral function 
of the w-meson has been discussed in detail in p6| . Here, the reaction lutt — > hit was found to cause a 
considerable broadening of the ui spectral function, leading to r w (Tc) ~ 7 r w (0). The corresponding 
photon spectral function is displayed in figure || (left panel) . 

At higher invariant masses 1 GeV < M < 2 GeV, irai annihilation is the dominant source of dileptons 
p7j , |28| . The vacuum vector and axialvector spectral functions become mixed to order T 2 with a strength 
2^/(6/^) (/tt ~ 93 MeV is the pion decay constant) due to their coupling to the pionic heat bath |2j| 
and should be degenerate at the point of chiral symmetry restoration. The effect of the a\ and higher 
resonances can then be approximately subsumed in a structureless continuum above 1 GeV |3^]. We 
practically implement the 7rai contribution by adding a flat 2-7T continuum to the resonance p meson 
spectral function that feeds into the photon spectral function, reminiscent of the perturbative plateau 
of qq annihilation. 

There is still considerable stopping of the interpenetrating nuclei at SPS energies, resulting in a net 
baryon density ps in the central rapidity region. At RHIC, measurements |3^| indicate that the proton 
over antiproton excess is small, implying that the baryons are distributed over a larger rapidity interval. 
Therefore, finite baryon density effects should not play such an important role at RHIC kinematics. 

For the evaluation of density effects which are relevant at SPS conditions, we use the results discussed 
in j3j| . There it was shown that in the linear density approximation, ny is related to the vector meson 
- nucleon scattering amplitude Tvn- 

Tl v (q°,q = 0;p) = H^ c - p B T VN (q), T VN {q) = - l - J d^x {N \T ^{x) f {Q)\N) , (17) 

with \N) being zero-momentum free nucleon states. In the following, we assume that the temperature- 
and density dependences of ny factorize, i.e. we replace n™ in eq.rtlT]) by the temperature-dependent 
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Figure 4: The photon spectral function R(u>) = (12ir/ui 2 ) Imfl(w) at finite temperature and ps = (left 
panel) and at T — and baryon density of normal nuclear matter, pb = pa = 0.17 fm -3 (right panel). 
For orientation, the qq line in the left panel shows the spectral function in the QGP phase with massless 
u- and d-quarks and m s = 150 MeV for s-quarks, neglecting a s -corrections. 



Hv(T) and leave Tvn unaffected. This amounts to neglecting contributions from matrix elements 
such as (ttN \Tj ll (x)j ll (0)\TTN) (nucleon-pion scatterings where the pion comes from the heat bath). 
Furthermore, this approximation does not take into consideration a possible T-dependent pion or nucleon 
mass. Some effective models suggest that, near the phase transition, the nucleon mass follows the 
behaviour of the chiral condensate (ipip) and drops abruptly as the quarks lose their constituent masses. 
Such modifications of particle properties may have a considerable impact on the spectral functions. 
However, since the temperature range over which the dropping takes place is narrow, we expect such 
effects not to leave distinct signals in the dilepton spectra which are only sensitive to the integrated time 
(and hence temperature) evolution of the system. 

The photon spectral function at finite density and zero temperature is depicted in figure |5| (right panel). 
The interaction with nucleons causes a strong broadening of the p meson down to the one pion threshold, 
leading to a complete dissolution of its quasiparticle peak structure. The modifications of the to and <fi 
meson spectral distributions are more moderate: The mass of the ui drops by about 100 MeV at normal 
nuclear matter density, and its width increases by a factor of about 5, whereas the <j> mass stays close to 
its vacuum value, accompanied by a ninefold increased width. 

To summarize, the most prominent changes of the photon spectral function, when compared to the 
vacuum case, arise from the broadening of the p due to finite baryon density effects and the broadening 
of the u> due to scattering off thermal pions. The (j) meson retains its distinct peak structure even 
under extreme conditions of density and temperature. Very close to Tq, however, these results based on 
perturbative calculations, are not expected to be reliable. 



4.3 After freeze-out contributions 

At the freeze-out stage, there are still vector mesons present. These will decay with their vacuum 
properties on their way to the detector and add to the dilepton yield from the previous thermalized phase. 
The invariant mass region below approximately 400 MeV is mainly filled by the Dalitz decays of the 
vector mesons. We take these contributions from the experimental analysis of the CERES collaboration 
for SPS conditions. Since the PHENIX acceptance starts only above 1 GeV, the Dalitz decays do not 
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play a significant role at RHIC. 



For the calculation of the direct decay of a vector meson V into a lepton pair we start with the following 
formula: 

Here, Tj and V/ are the fireball temperature and volume, respectively, at freeze-out. After Tf, the freeze- 
out time, all medium effects arc switched off, so the vacuum photon spectral function Ry(M,T = 0) 
determines the rate. The corresponding momentum distribution is given by the thermal Bose function, 
evaluated at the freeze-out temperature Tf. However, the absolute number N of vector mesons as a 
function of time is not a constant: since the mesons decay and there is no thermal recombination, N 
decreases exponentially like exp (— (r — T f) / ■ The vacuum life time of the vector meson V 

under consideration is denoted by Tq , and 7(9) accounts for time dilatation effects on particles with 
finite three-momentum: 

1 q° 

After the time integration we end up with: 

dN 1 



jT Q R{M, T = Q)V f J d 3 q f B (q°, T f ). 



dMd-q A?7l27r 

The averaged space-time four- volume that is available after freeze-out is therefore V/Tq , as anticipated. 
The integral over d 3 q now yields the freeze-out particle density n(M) — N(M)/V of the virtual photons 
as a function of invariant mass. Note that the information on the vector mesons remains entirely in 
the photon spectral function. With the factor Vf we obtain the total number of photons at freeze-out. 
When weighted with Ry{M), this gives the dN/dM distribution for the process meson — ► 7* — > e + e~. 

We have also checked that a commonly used Breit-Wigner ansatz of the form 

~dET ~ [Vfn(Tf > M)] ^ Bv (M 2 - m 2 v ) 2 + M 2 T V (M) 2 (18) 

yields almost identical results. Here, Ty(M) stands for the M-dependent decay width, £ is the normal- 
ization factor and By the branching ratio for the leptonic decay of the corresponding vector meson. 



4.4 Drell-Yan and charm contributions 



At invariant masses M > 1 GeV the Drell-Yan mechanism, i.e., hard quark-antiquark annihilation 
qq — > £ + £~ at leading order (LO), constitutes another source of (non-thermal) dileptons. Its differential 
LO cross section in a nucleus(Ai)-nucleus(yl2) collision reads 



da{A 1 A 2 ) %-ko 2 



9Ms ^ q 

q 



dydM 



(19) 



where </s denotes the cm. energy of the nucleon-nucleon collision and 3:1,2 = M/^fs exp(±y) is the 
momentum fraction of the beam and target parton respectively. 

The Drell-Yan cross section (|l|) is computed using the LO MRST parameterization [Q for the parton 
distributions /f(a;,^ 2 ) evaluated at the hard scale /j, 2 = M 2 . However, it has been checked that using 
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different LO sets (e.g., CTEQ5L |§ or GRV98LO [§6)) affects the resuits by only 10 % at SPS and 20% 
at RHIC energies. To account for higher order corrections, we multiply the LO expression (jll]) by a, K 
factor K — 2 fitted from p-p data |37j . Finally, nuclear effects like shadowing or quark energy loss are 
expected to suppress the Drell-Yan yield by about 30-50% |3^, [39|. Since these effects are still poorly 
known quantitatively, we neglect them and consider our Drell-Yan rate as an upper limit on the actual 
rate. 

Using the Glauber model of multiple independent collisions, the average dilepton multiplicity in a A1-A2 
collision at impact parameter b is given by 

dydM ( b )- T ^A 2 (b)xK dydM , (20) 

where the normalized thickness function Ta 1 a 2 Q ) ) is computed assuming the standard Woods-Saxon 
nuclear density profile. The Drell-Yan pair multiplicity (|2^) is then averaged for the 30% and 6% most 
central collisions to be compared with CERES and PHENIX data, respectively^. 

Another source of dileptons in the high invariant mass region consists of semileptonic decays of charmed 
mesons. Whereas earlier calculations found a considerable yield from open charm exceeding the thermal 
radiation Gfl], the subsequent inclusion of medium effects like energy loss led to a suppression of the 



dielectron rate and made it comparable to or even lower than the Drell-Yan yield Jl2|, 41 1. Since the 
Drell-Yan contribution plays only a subdominant role in the following, we will not explicitly include the 
open charm contributions. 



5 Fireball model 

5.1 General properties 

We do not aim at the detailed description of the heavy-ion collision or its subsequent expansion on an 
event by event basis or by a hydrodynamical simulation. Instead we use a model of an expanding fireball 
which enables us to test different scenarios with different parameter sets in a systematic way, so as to 
gain insight into the time evolution of the strongly interacting system. 

We assume that the physics of the fireball is the same inside each volume of given proper time r, thus 
averaging over spacial inhomogenities in density and temperature. The volume itself is taken to be an 
expanding cylinder, in which the volume elements move away from the center in order to generate the 
observed flow. There is no global Lorcntz frame in which thermodynamics can be applied. As the fireball 
expands, volume elements away from the center are moving with large velocities and are subject to time 
dilatation when seen in the center of mass frame of the collision. We assume a linear rise in rapidity 
when going from central volume elements to the fireball edge along the beam (z)-axis and the transverse 
axis. As the velocities along the z-axis are typically large (up to c) as compared to transverse motion 
(up to 0.55 c) for SPS and RHIC conditions, we make the simplifying assumption that the proper time 
is in a one-to-one correspondence to the z-position of a given volume element, thus neglecting the time 
dilatation caused by transverse motion. The whole system is assumed to be in local thermal (though 
not necessarily in chemical) equilibrium at all times. 

Given this overall framework, the volume expansion of the fireball is governed by the longitudinal growth 
speed v z and the transverse expansion speed vj_ at a given proper time. These quantities can be 

§)ln addition to these centrality cuts, we need to rely on further assumptions to take properly into account the acceptance 
of these experiments. Therefore we shall assume the generic form du/dp± oc pj_/(l + (p± /po) 2 ) 6 (po = 3 GeV) for the pj_ 
dependence of the DY process fjVm- Furthermore, neglecting corrections due to the intrinsic k± of the incoming partons, 
the angular distribution is taken to be dN/dQ oc 1 + cos 2 9 where 8 is the angle between the lepton and the beam axis. 
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determined at the freeze-out point and correspond to the observed amount of flow. However, flow is 
measured in the lab frame and needs to be translated into the growth of proper time volume. We use 
a detailed analysis of the freeze-out conditions for central Pb-Pb collisions at 160 AGeV Ej| to fix the 
endpoint of the evolution. The initial state is constrained using the overlap geometry of the colliding 
nuclei. The expansion between initial and freeze-out stages is then required to be in accordance with 
the EoS as determined from the quasiparticle model described in section [?]. The resulting model serves 



as the basic setup, its extension to different beam energies and collision centralities is discussed in 5.5 



5.2 Initial and Freeze-out conditions 

We use the data set bl of |Q as the endpoint of our fireball evolution. The data set has been obtained 
by a simultaneous fit of the fireball emission function to hadronic spectra and HBT radii, so as to 
disentangle the contributions from flow and temperature to these quantities. 

The fireball is characterized by a transverse box-shaped density distribution with a radius Rb = 12.1 fm. 
This corresponds to a root mean square radius of R( ms — 8.55 fm (= Rb/V%)- The average transverse 
expansion velocity is found to be 0.5 c and the temperature at freeze-out is Tf = 100 MeV. 

Looking at the longitudinal expansion, a velocity vt, « 0.9 c at the fireball front is needed in ±z direction 
in order to account for the observed shape of dN c h/dy distributions. 

The initial conditions in transverse direction can be fixed by the overlap geometry. Here, Rq ~ 4.5 fm 
and v± = for central collisions. The initial longitudinal size of the fireball is related to the amount 
of stopping of the matter in the collision and the time To necessary for the formation of a thermalized 
system. There is no direct information on these two quantities. However, pQCD calculations indicate 
To ~ 1 fm/c for SPS conditions and shorter times at RHIC. The initial longitudinal velocity v\ can be 
inferred from hydrodynamical calculations. Since the initial state in principle determines the final state 
(if the EoS is known) , one can fit this parameter to the observed dN / dy spectra. This procedure points 
at v\ w 0.5 c. The initial longitudinal size of the system at proper time To is then calculated by the 
intercept of z(t) = v\ ■ t with the t = y/t 2 — z 2 — tq line. 



5.3 Volume evolution 

Using the available information on initial and freeze-out conditions and the EoS of the system, we are 
now able to reconstruct the evolution of the fireball volume in proper time: 

(1) The EoS translates into a temperature (and hence t) dependent acceleration profile a = const. ■ 
p(T)/e(T) that exhibits a soft transition point at T — Tq- This can be understood as follows: the 
thermal vacuum energy B(T) adds to the energy density e but is subtracted from the pressure p. As 
evident from figure I B(T) gets large near the phase transition and correspondingly the ratio p/e 
drops. In the interpretation of B(T) as a 'bag pressure', the pressure of the QGP becomes partly 
counterbalanced by the confining bag pressure. Comparing initial and final conditions, we need both 
longitudinal and transverse acceleration to account for the velocity gain, and we keep the possibility of 
having two different constants cj_ and c z which relate transverse and longitudinal accelerations to p/e. 
In practice, the temperature evolution with r is calculated starting with a trial ansatz and iterating the 
result to obtain a self-consistent solution. 

(2) Starting with an ansatz for the radial expansion velocity, 

«x(r) = jTdr'cx^ (21) 
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and 

R(r) = R + f T f T dT'dr"c ± ^l (22) 
Jo Jo e\ T ) 

where Rq is the initial overlap rms radius of the collision region, we fix the two unknown parameters c± 
and Tf by the requirement that R(tA = R( ms and v^(ry) = Vj_. 



(3) For the longitudinal expansion we follow the motion of the fireball front in the center of mass frame 
and use the expressions 

v z (t)=vi+j t dt'c z ^l (23) 

and 

z(t) = z + vl ■ t + /* / dt'dt"c z ^Q-. (24) 
Jo Jo e l t ) 

The parameters c z and tf (freeze-out in cm. frame) are fixed by v z (tf) — v( and z(tf) to lie on the 
t = Tf line, with r/ determined from the radial expansion. 



(4) The movement of the fireball front in the cm. frame can now be translated into the growth of 
volumes in proper time by finding the intercept of z(t) and r = const, and calculating the pathlength 
along this curve of fixed proper time. 

In order to construct a self consistent evolution model, we need to find the proper ratio p(T '(r)) / 'e(T '(r)) 
corresponding to the volume expansion. The evolution of T(r) is dealt within the next section. For 
the QGP phase, the ratio is then determined within the quasiparticle model. Unfortunately, no reliable 
information on p/e is available for the hadronic phase. We can, however, deduce the value of this ratio 
at both Tq and Tf, and interpolate between these limits. As the lattice indicates neither a sharp drop 
in e nor in p when approaching Tq from above, it appears reasonable to assume that the p/e ratio 
stays small in the vicinity of the phase transition even in the hadronic phase. On the other hand, 
the decoupling of the system at freeze-out implies that interactions between its constituents become 
unimportant, therefore we recover the ideal gas limit at T — > Tf, where standard thermodynamics 
predicts p/e — 1/3 (for massless particles). We now interpolate linearly between these two values to 
cover the temperature region inbetween. 



5.4 Temperature evolution 

The temperature profile T(t) of the fireball is uniquely determined by the assumed condition of isentropic 
expansion once the volume expansion is known. In order to derive the temperature T at a given time, 
we calculate the entropy density via 

s(t) = S /V(r), (25) 

where So is the total entropy of the fireball. The relation between the entropy density s and the fireball 
temperature T is fixed by the EoS. By inverting this (unique) relation numerically, we finally obtain the 
temperature profile 

We infer the EoS of the QGP phase from the thermal quasiparticle model |l!| adapted to lattice QCD 
results. Unfortunately, determining the behaviour of the hadronic phase is not quite as easy for a num- 
ber of reasons. First, lattice results usually employ quark masses which are far too large, leading to 
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pion masses m l ^ tt > 3 mP hys . The resulting thermal suppression of these degrees of freedom causes a 
considerable discrepancy of the lattice EoS with respect to that of the interacting pion gas or the free 
hadronic resonance gas. On the other hand, the description in terms of a non-interacting or perturba- 
tively interacting system is bound to fail near the phase transition. Only in the very final stages, near 
the freeze-out, interactions cease to be important and one may assume that a non-interacting system 
describes the situation adequately. 

We parametrize our insufficient knowledge close to Tc by interpolating smoothly between two regimes, 
the QGP quasiparticle result for T > Tc and the noninteracting hadronic resonance gas result for 
T < Tf. This approach is supported by the general shape of the EoS on the lattice for two light quarks 
and one heavy quark, which resembles a weak first-order or second order phase transition, or even a 
smooth crossover 
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Figure 5: The temperature dependence of the entropy density s(T) in the RHIC kinematical scenario 
as compared to the ideal quark-gluon gas and the massless noninteracting pion gas limit (dotted line). 
The three relevant regions used in the model calculation are given as ideal hadronic resonance gas (solid 
line), interpolation (dash-dotted line) and thermal quasiparticle model (dashed line). 

The situation is shown for the RHIC scenario in figure [5| and compared to the ideal quark-gluon gas 
and the hadronic resonance gas approach used in other works (e.g. Q|). One can clearly observe 
that the deviations in both the quark-gluon phase and the hadronic phase from the ideal gas are large, 
amounting to a factor of more than two for the entropy density near Tq. Keeping in mind that a smaller 
entropy density translates into a higher temperature in an isentropic expansion, we conclude that our 
model predicts a prolonged lifetime of the QGP phase as compared to the ideal gas ansatz, whereas the 
lifetime of the hadronic phase is reduced somewhat. 

At SPS, the measured ratio of protons over antiprotons indicates a partial stopping of nucleons during 
the collision phase, resulting in an excess of quarks over antiquarks in the fireball region jl4j. Since 
baryon number is conserved, this implies the existence of a T-dependent baryochemical potential, which 
in turn translates into a fugacity factor As = exp[^-] multiplying the entropy density created by baryons. 
Furthermore, the thermal yield of pions in the fireball is not enough to account for the observed number 
of pions. This can be compensated by the introduction of a chemical potential for pions (and kaons), 
which in turn influences the EoS. 

By imposing entropy and baryon number conservation, the evolution of /is can be followed through 
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the fireball expansion. In practice, the evaluation at each T rests on the assumption of the system 
being an ideal hadronic gas, which we believe is unreasonable near the phase transition. Fortunately, 
Tf ~ 100 MeV is far enough distant from the phase transition, and that is where we fix the entropy. 
Furthermore, we expect pions to be the dominant thermally active degrees of freedom. We use a pion 
chemical potential /j, n (Tf) = 123 MeV as obtained in p5j, which gives the correct total pion multiplicity 
when evaluated with the fireball freeze-out volume. 

Correcting the contributions from the different particle species for the corresponding fugacity factors, 
we obtain a point of the s(T) curve where our interpolation starts. This situation is shown in figure ^. 
One observes that, unlike the RHIC scenario, the entropy density under SPS conditions is larger in the 
range Tf < T < Tc, resulting in a faster dropping of temperature during the hadronic expansion phase. 
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Figure 6: . Temperature dependence of the entropy density for SPS conditions, including the interpola- 
tion to the fugacity corrected value at freeze-out. As the entropy density below Tf is irrelevant for the 
fireball evolution, the interpolation is stopped at this point. 

Having now specified the relevant part of the EoS of the system, the temperature evolution of the fireball 
is uniquely determined by the total entropy So and the evolution of the volume V(t) with proper time. 
The total entropy can be obtained by measuring charged particle multiplicities N + and N~ in suitable 
rapidity bins and calculating 

N+ - N~ 

D * = iv+T^F' (26) 

The quantity Dq stands for the inverse of the specific entropy per net baryon S/B, and the product 
Dq(S/B) roughly measures the entropy per pion For SPS collisions at 160 AGeV, we find S/B = 26 
for central collisons. For RHIC 6% central Au-Au collisions at 130 AGeV, the specific entropy S/B = 130 
is substantially higher due to the larger particle multiplicity and the smaller net baryon content in the 
central region. For beam energies of 200 AGeV that are of interest here, not enough information on 
dN/dy spectra is available at this moment, so we utilize the predictions from a thermal model calculation 
[ lira . With the ratios p/p — 0.75, p/ir~ — 0.09 and K~ /ir~ = 0.15, we obtain the specific entropy as 
S/B ps 250, as already estimated (albeit with different parameters) in p3fl. 



17 



5.5 Variations with beam energies and centralities 



It is mainly the detailed information on the final and initial state which enables us to reconstruct the 
fireball evolution in the case of the 160 AGeV central collision scenario. Unfortunately, no such detailed 
analysis is available for the 30% central PbAu collision scenario for which dilepton data have been taken. 
The same holds for the 40 AGeV dilepton data. It is therefore necessary to extend the framework 
established so far to different centralities and beam energies by a suitable re-scaling of key quantities 
characteristic of the evolution. 

First of all, the total entropy deposition in the fireball region must be different when going to more 
peripheral collisions or lower beam energies. We assume that the entropy per baryon scales with the 
number of negatively charged hadrons observed in the final state and take the SPS value of s/ps = 26 
as a baseline. When going to peripheral collisions, we assume that this value is still a good estimate. 
Here, the total entropy in the system is naturally reduced because the number of nucleons participating 
in the collision is smaller. 

Second, the overlap geometry is different in peripheral collisions, resulting in a smaller initial fireball 
volume. Here, we neglect details of the transverse geometry and keep parametrizing the fireball volume 
as a cylinder. Its initial transverse area is adjusted to the value of the calculated overlap area, hence the 
cylinder radius is reduced as compared to central collisions. 

For the 40 AGeV data, we take the total entropy to be about half the entropy at 160 AGeV because 
N 160 as 600 and iV 4 ° sw 310. The initial energy density, estimated by Bjorken's formula from dN/dy 
at midrapidity, is about 2/3 of the value at 160 AGeV, but still well above the critical energy density 
required to form a QGP plasma. Looking at the final state of the fireball, HBT analyses of NA49 JE^] 
indicate that the radius parameters are very similar at 40 and 160 AGeV, suggesting that the reaction 
dynamics do not significantly change in this energy region. 

Changes in the initial entropy deposition result in a correspondingly different fireball evolution and in 
general a different freeze-out state. Freeze-out occurs when the mean free path of particles exceeds the 
dimensions of the fireball. As pions are the most abundant particle species in the fireball, we assume that 
the pion density determines the mean free path A of particles in the medium. The freeze-out condition 
reads therefore 

a\p„ = 1, (27) 

where a stands for a typical hadronic cross section. As the observed HBT radii for 40 and 160 AGeV 
appear rather similar, we do not expect the total freeze-out volume to change more than a factor two. 
This is thermodynamically consistent with fixing the freeze-out temperature Tf — 100 MeV for all SPS 
scenarios. As A scales ~ \/V, this is a sensible ansatz — the resulting freeze-out pion densities are 
very similar when looking at (p7|). This is still true if we take modificatios of the pion density by a 
pion chemical potential p w into account — about the same value of \i v is needed in all SPS scenarios to 



account for the observed total pion yield (see section 5.7) 



The situation is different at RHIC. Here, no large pion chemical potential appears to be necessary to 
account for the observed total yield and therefore the thermal pion density at a given temperature is 
lower, leading to a higher freeze-out temperature with (p7|) . We find that an expansion scenario with Tf 
= 130 MeV fits both the observed total particle yield and leads to the correct freeze-out pion density. The 
freeze-out volume has now to be adjusted accordingly. In order to do this consistently, we simultaneously 
modify both the final state flow velocities and the fireball radius as compared to our standard scenario. 
As these quantities result to a first approximation from an accelerated motion, reducing the radial flow 
by a factor / implies a reduction of the longitudinal flow difference between initial and final state by 
the same factor / and a reduction of the freeze-out radius by f 2 . Note that a combination of decreasing 
geometrical radius and flow velocity tends to keep the HBT radius constant. 
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The fireball evolution is now calculated as described above, using re-scaled total entropy, initial radius, 
final state longitudinal and transverse flow and freeze-out radius as new inputs. The parameter sets 
obtained for the different fireball scenarios are summarized in the appendix. 

5.6 Results and discussion 

The resulting volume evolutions for 40 and 160 AGeV are plotted in figure [?]. Note that these curves 
correspond to the volume at proper time r, which is larger than the geometrical size of the fireball in 
the cm. or the lab frame. The temperature evolution is shown in figure ^. We observe that the QGP 
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Figure 7: The volume expansion for SPS conditions. Shown are the curves for central 160 AGeV 
collisions (solid line) and central 40 AGeV collisions (dashed line). 

phase lasts about half of the total fireball lifetime for 160 AGeV, and considerably less for 40 AGeV. 
As the fireball evolution resulting from our model differs somewhat from the results obtained by other 
groups, we feel that some clarifying remarks are in order. 

In our quasiparticle model there is no mixed phase of coexisting hadronic gas and QGP. Based on the 
observation that there is no strong first order transition visible in the lattice data, the EoS right and left 
of Tc match smoothly. Inserting the EoS of an ideal quark-gluon gas instead, a large gap in the s(T) 
diagram is created. The resulting latent heat (AS)Tc generates a mixed phase of considerable duration 
(5 - 10 fm), as found in previous approaches. 

By construction of the model, hadronic observables derived from the expanding fireball, e.g. rapidity 
distributions of charged particles, p t -spectra, HBT radii and observed particle numbers, are described 
adequately. This is an important constraint of the model. There is very little freedom left to tune the 
fireball evolution. Arguably, the EoS inferred from the quasiparticle model is the weakest point in the 
chain of arguments leading to the complete picture. However, bulk quantities such as freeze-out proper 
time and freeze-out volume are hardly influenced by the EoS or the phase transition temperature. They 
can be found even by assuming constant acceleration at all times. 

The initial temperatures are quite large (~ 300 MeV for SPS central conditions, 420 MeV for RHIC) 
and are uniquely determined by the total entropy, the initial volume and our EoS for the quark-gluon 
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Figure 8: Time evolution of the temperature for SPS and RHIC conditions as obtained with the fireball 
model described in this work. The critical and freeze-out temperatures are indicated by the horizontal 
lines. 



phase. To some extent these large values related to our choice of the formation time To = 1.0 fm/c (SPS) 
and To = 0.6 fm/c (RHIC) of the thermal system; if To is increased, the initial temperature decreases 
correspondingly. However, unlike the results found in 0], we do not observe any strong sensitivity of 
the dilepton yields to the choice of To: as the fireball expands, differences in the initial volume become 
increasingly unimportant. If the total evolution time of the QGP phase is small, these changes may 
well matter, but as the lifetime of the QGP phase in our model is comparatively large, the possible 
difference affects only a small fraction of the total lifetime in a region where the fireball volume (and, 
correspondingly the dilepton yield) is small anyway. 

Comparing with the Bjorken estimate of the initial energy density, one should keep in mind that the 
extrapolation from the final to the initial state is different in the present approach. In the Bjorken 
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scenario, no longitudinal acceleration is present, therefore the mapping of final state rapidity distributions 
to initial state spatial distributions finds a larger initial volume than our scenario. If we assume no or 
only small longitudinal acceleration in order to compare the two approaches, we find initial temperatures 
between 220 and 240 MeV, consistent with the Bjorken estimate. As this assumption is incompatible 
with a freeze-out at 100 MeV and the measured freeze-out geometry, we disregard it. 

We observe a prolonged lifetime of the QGP evolution phase of the fireball as compared to the results 
obtained by other groups. Recall that the apparent lifetime of the QGP, as observed in the lab frame, is 
larger than its lifetime in proper time r. This is a consequence of our volume evolution and the use of 
the more realistic EoS of our quasiparticle model, as opposed to that of the ideal quark-gluon gas. Near 
the phase transition, the entropy density in the quasiparticle model is about a factor two smaller than 
the one of the ideal quark-gluon gas, and it takes a correspondingly larger volume (and larger evolution 
time) to reach this entropy density in an isentropic expansion. 

The slower cooling time in the hadronic phase of the RHIC scenario as compared to SPS conditions can 
in essence be traced back to the fact that the entropy density under SPS conditions is enhanced by the 
fugacity factors whereas this is not such an important effect for RHIC. 

5.7 Chemical composition 

As already mentioned, the total thermal pion yield of the fireball at freeze-out is not enough to account 
for the observed number of pions, and a large chemical potential of p n = 123 MeV has to be used to com- 
pensate ijlsl . The common point of view is that particle numbers are fixed at a chemical freeze-out point 
and the absence of inelastic reactions introduces then a chemical potential. The statistical model is 
enormously succesful in describing the finally observed ratio as a fit of chemical freeze-out temperature 
and baryon density, using the free properties of particles. However, in the hot and dense environment 
of a fireball, in medium modifications of particle masses and widths are likely to occur, changing the 
amount of particle production at the freeze-out temperature. This is not in direct contradiction to the 
apparent success of the statistical model, as a different set of T and fx might be able to explain the 
observed ratios once the modifiactions are taken into account. Furthermore, the absence of all inelastic 
reactions, leading to decay processes only, is most likely an oversimplification for some particles, e.g. for 
the p <-> 7T7T system. If one knew the particle properties and the particle abundancies at hadronization 
reliably, the correct way to proceed would be to set up a system of coupled rate equations and follow 
the various decay and recombination processes to the thermal freeze-out. Unfortunately, theoretically 
there is no way to assess these properties reliably near the phase transition, as this is a nonperturbative 
problem. Keeping these issues in mind, we will therefore parametrize the hadrochemistry of the fireball 
on a phcnomcnological basis instead of aiming for a detailed solution of the problem. 

We adopt the following approach: A pion chemical potential is introduced ad hoc and fitted at the 
thermal freeze-out point (p, n = 123 MeV |45[]) to the total pion multiplicity. It is assumed to decrease 
with temperature linearly up to Tq where it vanishes, corresponding to a situation where resonance 
decays continuously feed pions into the system. All other necessary chemical potentials are adjusted 
such as to reproduce the observed pion to particle ratios at all times. 

6 Dilepton invariant mass spectra 

Once the time evolution of the fireball is given in terms of the temperature T(r), the baryon density 
p(r) and the volume V(r), and with the knowledge of the photon spectral function, we have all the 
necessary ingredients to calculate dilepton rates using eq. (||). We fold the result with the acceptance 
of the CERES and the PHENIX detectors, respectively, and average over the rapidity region covered by 
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these two experiments. The so-called 'hadronic cocktail', dileptons produced after freeze-out through 
various decay processes, with the exception of vector-meson decays, is then added. This contribution fills 
the region of very low invariant masses (M < 150 MeV). The dilepton yields resulting from direct vector 
meson decays after freeze-out, as described in subsection 4.3, and the Drell-Yan yield from subsection 4.4 
are added to the hadronic cocktail, taking into account the limited kinematic acceptance and resolution 
of the detector. 



6.1 SPS data at 40 and 160 AGeV 

We start with a discussion of the SPS CERES/NA45 experiment, treating 40 AGeV and 160 AGeV data 
in parallel. The final results for the dilepton invariant mass spectra are shown in figs. |^ and [IT^. 

Our calculation reproduces the overall spectrum of the 160 AGeV CERES data quite well. It overesti- 
mates the rates somewhat around invariant masses of 200 to 400 MeV and achieves a good description 
in the region above 400 MeV up to 1.8 GeV. Recall that our QGP model rate constitutes only a lower 
limit on the actual rate because it neglects the radiation from non-partonic (cluster) degrees of freedom 
above the critical temperature. Bearing in mind that the region above 1 GeV is mainly populated by 
dileptons originating from the QGP phase, as evident from the left panel of figure |^, there might still 
be additional radiation close above Tq arising from hadronic clusters embedded in a QGP environment. 
The Drell-Yan contribution is non-negligible, but still outshined by the QGP by a factor of 3. 

Changes in the spectra of the vector mesons indicate tendencies towards chiral symmetry restoration, 
so the right panel of figure [)] shows the contributions of p, u> and <j) mesons separately, not including 
the after freeze-out yield. The p meson loses its quasiparticle structure entirely due to strong collision 
broadening at finite density, and fills the whole low mass region between the two pion threshold and 
~ 800 MeV. The to meson, a sharp peak in the vacuum, broadens at finite temperature due to the 
thermal scattering process lot: — * tttt . Furthermore, the mass shift at finite baryon density smears the 
remaining peak structure considerably. Effectively, the remaining signal arises from the direct decays of 
u mesons after freeze-out. The <f> has become a spread-out but still visible resonance structure, showing 
only moderate broadening at finite temperature and baryon density. It might therefore be a suitable 
candidate for gauging the strength of vector meson modification. 

To test the modelling of the vector meson spectra, we calculate the total number of u) and (f> meson, 
suitably averaged over their medium-induced spread in invariant mass, as 

0.65 GeV 

and 

1.1 GeV 

, , 1 f CpN^e , 

^ = N7 h J dM ^kr - L8 ' 10 

0.95 GeV 

Comparing with numbers from a statistical model calculation, (NJ) = 4- 10~ 7 and (N^) = 2.2-10~ 7 [j54| , 
we indeed find reasonable agreement. Note that the relatively large u meson yield is primarily caused 
by the pion fugacity factor [exp(/i 7r /T/)] 3 at freeze-out, which reflects the enhanced feeding through the 
37r — * oj process. 

Going from 160 AGeV to 40 AGeV beam energy, analyses of HBT radii and transverse radial flow ]5^| 
indicate that the reaction dynamics do not change dramatically, therefore we do not expect drastic 
differences in the rate. Indeed, the data at 40 AGeV look similar to the 160 AGeV case, and the 
calculated rate, shown in figure ITol, also bears this similarity and achieves a good fit without tuning 



22 




Figure 9: Left: Dilepton invariant mass spectra (see eq. (0)), normalized to dN^/dr/ = 250, in units 
of (100 MeV)"\ for the SPS CERES/NA45 Pb(158 AGeV)+Au experiment §§. Shown are the data, 
the total rate, the cocktail contribution including the after freeze-out decays of vector mesons, the QGP 
contribution and the Drell-Yan yield. Right: Contributions from p-, uo- and 0-mesons (excluding after 
freeze-out yield) shown separately, assuming perfect detector resolution. 



Pb(40 AGeV) + Au 




invariant mass M [GeV] invariant mass M [GeV] 



Figure 10: Left: Dilepton invariant mass spectra (see eq. (|3|)), normalized to dN c h/drj — 210, in units of 
(100 MeV)" 1 , for the SPS CERES/NA45 Pb(40 AGeV)+Au experiment [g. Shown are the preliminary 
data, the total rate, the cocktail contribution including the after freeze-out decays of vector mesons and 
the QGP contribution. Right: Contributions from p-, u>- and 0-mesons (excluding after freeze-out yield) 
shown separately, assuming perfect detector resolution. 



the setup of the model. Since the initial temperature is lower and the QGP phase shorter in the 40 
AGeV case, the partonic dilepton contribution is obviously much smaller, but nevertheless still present. 
Owing to the greater initial baryonic density, the in-medium modifications of the vector mesons become 
more pronounced, most prominently visible in the uj meson channel. Its downward mass shift drags the 
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peak structure along the time evolution of the fireball, creating a small bump on top of the completely 
dissolved p meson that fills up the low-mass region again. Its yield after freeze-out constitutes a visible 
signal that may be experimentally observable with suitable energy resolution. The tf> meson contribution 
clearly sticks out above the smooth p meson 'continuum'. To conclude, we find no distinct differences 
in our calculation for the two beam energies probing dilepton production at SPS so far, in accord with 
experimental findings. This indicates that the general setup of our model is fairly robust. Future data 
at 20 and 80 AGeV will aid to test this statement. 




Figure 11: Left: Dilepton invariant mass spectra for transverse momenta of the e + e pair p t < 500 MeV 
for the SPS CERES/NA45 Pb(158 AGeV)+Au experiment §5|. Shown are the data, the total rate and 
the cocktail contribution. Right: Same for p t > 500 MeV. 




Figure 12: Left: Dilepton invariant mass spectra for transverse momenta of the e + e pair p t < 500 
MeV for the SPS CERES/NA45 Pb(40 AGeV)+Au experiment §7|. Shown are the data, the total rate 
and the cocktail contribution. Right: Same for p t > 500 MeV. 

The fact that we moderately overestimate the data in the region between 200 and 300 MeV invariant mass 
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requires a comment. Since this range is dominated by the low- mass behaviour of the p meson spectral 
function at finite density, this behaviour may indicate that the influence of finite three-momentum on 
the spectra in that very region is non-negligible. Consider the f 60 AGeV and 40 AGeV data taken for 
different transverse momenta pt < 500 MeV and pt > 500 MeV, shown in figures [Tl] and [l2|. We observe 
that the general shape of the data pattern is well described by the calculation for both p t regions and 
both beam energies. However, for the high p t case, the calculation again overshoots the data in the 
low mass region whereas this is not so in the low p* case. This can be traced to the use of the spectral 
functions for three-momentum equal to zero in order to describe the photon spectral function in the 
hadronic phase. This approximation is of limited validity at high p t where the spectral function ought 
to become smaller than in our approximation. Therefore we expect improved agreement with the data 
for the low invariant mass region both at 40 and f 60 AGeV once this effect is taken into account properly. 

Our setup of the fireball model enables us to gain additional detailed insight into the time evolution of 
the dilepton yield. This is demonstrated in figure Here the different stages of the fireball evolution 
leave distinct marks in the time-resolved dilepton yield. For early times, dileptons come entirely from 
the qq quasiparticle annihilation processes. The movement of the invariant mass threshold reflects 
the temperature dependence of the quasiparticle mass which decreases near the phase transition at 
r ~ 7 fm/c. One observes that, in spite of the growing fireball volume, the contributions from later 
timeslices to the total yield become progressively less important until the hadronic phase takes over at 
r > 7 fm/c. This surprising behaviour is enforced by the confinement factor C'(T) which reduces the 
thermodynamically active degrees of freedom significantly near the phase boundary. The system then 
enters the hadronic evolution phase without going through a mixed phase. The most prominent feature 
is the rapid filling of the low invariant mass region through the density-broadened p meson spectrum 
which ends up as an enhanced pion continuum. The <f> meson starts contributing its characteristic peak 
and, as the system cools down further, the u> meson begins to emerge, albeit as a broad structure. Note 
that while the hadronic contributions fill primarily the low invariant mass region, their yield above 1 
GeV is negligible in this late evolution phase. 



6.2 RHIC at ^ = 200 AGeV 

For the RHIC scenario, thermally generated dileptons are dominant. Measurements of proton ratios at 
\/s = 130 AGeV indicate that the central collision region remains almost net-baryon free, compared to 
SPS energies. Within a statistical thermal model, the particle ratios are accounted for by a small baryon 
chemical potential of about 50 MeV at chemical freezeout [Q . At 200 AGeV, this value is predicted to 
be even smaller. Effects of finite baryon density are therefore almost absent and consequently both the 
p and the 4> are expected to show up in the spectrum as pronounced structures, whereas there should be 
no clear trace of the in-medium to due to its strong thermal broadening. Contributions from Drell-Yan 
processes, which dominate in the very high invariant mass region, are an order of magnitude smaller. 

The corresponding prediction for the dilepton yield at 200 AGeV, including the schematic acceptance 
for the PHENIX detector, is shown in figure ||. The lu and <f> meson resonances clearly stick out over the 
smooth p meson and QGP contributions. Although PHENIX will only start to measure at M > 1 GeV, 
it may be possible to resolve the <f> meson peak. However, a significant part of the peak strength is built 
up by the after freeze-out contributions, making it difficult to disentangle the in-medium modifications 
on the hadrons. For M > 1.3 GeV, the dilepton spectrum is dominated by thermal QGP radiation, 
outshining the hard Drell-Yan dilepton yield. 

Comparing our prediction for PHENIX with the one shown in rcf. fl43| , we find rough agreement of 
the rate for the low mass region below ~ 1 GeV. Although the dilepton yield from the QGP phase is 



suppressed in our case by the factor C(T) 2 (cf. section 4.1), we still find an enhancement of a factor 
of about 4 in the range 1.3 — 2.5 GeV over the rate in |43] that employs a (perturbative) chemical 
undersaturation model in the QGP phase. Owing to the non-perturbative nature of the QGP close to 
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M [GeV] 

Figure 13: Time evolution of the (integrated) dilepton yield for the SPS scenario. The vertical scale 
is the same as in fig. ||. The cocktail contribution has not been included in this picture, and perfect 
resolution was assumed. 

Tc, this model may not be applicable in that very region. Note that our rate also constitutes only a 
lower limit there, so that the actual rate may even be larger. High precision data will allow to pin down 
one or the other model. 
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Figure 14: Left: Dilepton invariant mass spectra, normalized to dNch/drj = 650 in units of GeV -1 , 
for the RHIC experiment PHENIX at ^/s = 200 AGeV. Shown are the total rate, the cocktail consisting 
of the after freeze-out decays of u) and <j> vector mesons, the QGP contribution and the Drell-Yan yield. 
Right: Contributions from p-, uj- and c/>-mesons (excluding after freeze-out yield) shown separately. 



6.3 Sensitivity to model parameters 

We would like to stress that the gross features of our model are set, once the parametrization of the 
fireball evolution has been matched to the hadronic observables and the EoS of both phases has been 
constructed in accordance with lattice QCD and empirical constraints. The remaining uncertainties, 
mainly about the initial state of the fireball, the thermal masses of the quasiparticles and the detailed 
shape of the EoS, do not alter the results substantially; they lead to only moderate or even weak 
dependence on those parameters. Fine-tuning is still possible, but only within the limits that retain 
consistency with the overall framework. 

We have investigated the sensitivity of the model to parameter changes in some detail for the SPS 
scenario at 160 AGeV. In order to get a theoretical 'error band', we investigated the extremes of our 
parameter ranges, a combination of parameters that yields the largest and the smallest possible QGP 
contribution. The resulting range is shown in figure [15| as a grey band, together with the data points 
and the curve from the previous section. 

Regarding the initial conditions, the largest uncertainty comes from the initial fireball formation length 
zq or, equivalently, the time tq. We let it vary from 0.5 fm/c to 2 fm/c, i.e. from fast to slow equilibration. 
Consider next the initial longitudinal flow velocity v\. A variation of this quantity from 0.3 c to 0.8 
c implies a modification of the pressure-driven accelerated motion in order to arrive at the same final 
velocity of 0.9 c, as determined by the rapidity distributions of the observed hadrons. Strong deviations 
from v\ ~ 0.5 c, however, lead to inconsistent values for the thermal freeze-out temperature TV. Finally, 
modifications of the critical temperature Tc influence mainly the relative weight of the contributions 
from the QGP phase and the hadronic gas phase to the dilepton yield and hence change the shape of the 
dilepton spectra. Lattice data on Tc suggest a range from 140 MeV for three massless, thermally active 
flavours to about 185 MeV for two massless flavours and a realistic pion mass. Due to the self-consistent 
modelling of the temperature and volume evolution, there is no simple one-to-one correspondence of 
these parameters to the dilepton yield. 

The upper limit of the grey band in figure |q now corresponds to the scenario with a small QGP 
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Figure 15: Dilepton invariant mass spectra for the SPS CERES/NA45 experiment at 160 AGeV. 
Shown are data (symbols), the total rate introduced in section |6.l| (solid line) and a band (shaded) that 
represents the range in the dilepton yield when varying the parameters To, v\ and T c . See text for details. 



contribution, i.e. small To, large v\ and high Tq- On the other hand, the lower limit includes a large QGP 
contribution with large To, small v\ and low T c . It is instructive to note that the shape of the spectrum 
changes only moderately above 1 GeV invariant mass within these extreme parameter variations. The 
first scenario, however, tends to overestimate the data in the region of the p peak, whereas the second 
scenario does not leave enough time for the hadronic phase to build up the e + e~ excess in the low-mass 
region between 200 and 800 MeV, effectively ruling out a large QGP contribution. 



7 Concluding Remarks 

We have calculated dilepton radiation from an expanding fireball created in ultra-relativistic heavy-ion 
collisions over a wide range of beam energies, from SPS 40 AGeV via SPS 160 AGeV to RHIC y^i = 
200 AGeV. 

We have explored the evolution of a fireball through the quark-gluon and hadronic phases as they are 
assumed to be characteristic of the QCD equation of state. Furthermore, we have emphasized the 
relevance of using input and constraints from lattice QCD thermodynamics in the analysis of dilepton 
signals from the expanding fireball. The fact that the measured dilepton rates could be reproduced using 
a fireball model fixed by hadronic information only is a nontrivial confirmation of its validity. Further 
data for different energies and/or nuclear systems will help to strengthen or weaken this position. 

High resolution measurements of the invariant mass spectrum around the region of the to and the 
</> mass can reveal information on the average density and temperature in the hadronic phase. The 
broadening and mass shift of the u> meson with temperature can serve as a 'thermometer' once the cocktail 
contribution is reliably assessed. A visible enhancement in the u> region beyond the rate originating from 
w-decays after freeze-out would hint at a temperature which is on average lower than assumed in this 
model and therefore point to a lower freeze-out temperature or alternative in- medium effects on hadrons. 
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On the other hand, the (f> meson signal, being nearly unmodified by temperature effects over a comparably 
long timespan, may be used as a 'standard candle' at RHIC conditions. 

Detailed knowledge of the slope of the invariant mass spectrum above 1 GeV gives an indication of 
the average temperature in the quark-gluon phase, serving as a constraint for the fireball dynamics in 
the QGP region. Its total normalization is useful to estimate the effects of a possible underpopulation 
of the quark and gluon phase space in early stages of the evolution, or of phenomenological models of 
the QGP, such as the confinement model presented in section |3|. It is remarkable that this confinement 
model is able to reproduce the data despite the suppression of dilepton radiation near the phase transition 
temperature. 

Measurements using different nuclear systems or different impact parameter can provide different relative 
weights of QGP and hadronic dilepton radiation to the observed yield. One may thus hope to disentangle 
these contributions and test the model assumptions above and below the phase transition separately, at 
least qualitatively. 

The present scenario supports the hypothesis that the quark-gluon phase is actually reached, at a 
transient stage, in heavy-ion collisions at CERN/SPS and (more pronouncedly) at RHIC. Within our 
outlined approach, a purely hadronic framework would not be successful in comparison with existing 
data. The reason is that the relative strength of hadronic contributions to the high mass region and 
the low mass region is quite different as compared to those of the QGP (cf . figure ^ . Purely hadronic 
scenarios which are able to account for the low mass region would necessarily fail in the high mass region 
and vice versa. Furthermore, fireball thermodynamics as described in this work leads to temperatures 
well above Tq for a broad range of initial conditions, though not necessarily to significant contributions 
to the dilepton yield from the high-temperature region. The upcoming measurements at RHIC can be 
expected to provide further interesting insights. 
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A Fireball properties 

In the following, we summarize key properties of the different fireball scenarios discussed in the paper. 
We characterize each scenario by the proper time for fireball formation, To, the phase transition time t c 
and the freeze-out time Tf. Note that the corresponding times in the cm. frame can be considerably 
larger, especially for RHIC conditions where volume elements are travelling near the speed of light at 
the fireball edge. Additionally, we quote the volume at freeze-out and the flow velocities at rms radius 
and maximum longitudinal extention v[. In order to compare the thermodynamic conditions, we 
provide s/p, the entropy per baryon. 



Note that the expansion is modelled in proper time. The lifetime of the fireball in the RHIC cm. frame, 
for example, is more than four times larger due to time dilatation at the fireball edge. 
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t [fm/c] 


t c [fm/c] 


T f [fm/c] 


Vf [fm 3 ] 




v z 


v f 


s/PB 


SPS 40 AGeV 


1.5 


4.0 


15.0 


7040 


0.36 


0.45 


0.75 


13 


SPS 160 AGcV 


1.0 


6.5 


16.0 


14.400 


0.53 


0.45 


0.9 


26 


RHIC 200 AGeV 


0.6 


10.0 


18.0 


99846.1 


0.56 


0.9 


0.9985 


250 



Table 1: Key parameters of the fireball evolution scenarios 
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